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The spin current and spin conductivity is computed through thermally driven stochastic process. 

By evaluating the Kramers equation and with the help of k.p method we have studied the spin 
Hall scenario. Due to the thermal assistance, the Kane model parameters get modified, which 
consequently modulate the spin orbit coupling(SOC). This modihed SOC causes the spin current 
to change in a finite amount. 


I. INTRODUCTION 

Today’s condensed matter research mainly relies on the study of the spin related issues of different materials. 
This put forward the concept of ” Spintronics” [1, 2], which unveils the importance of the spin degrees of freedom of 
electron for improved spin based devices. In this context, the development in the arena of semiconductor spintronics 
attracts the attention of many theoreticians as well as experimentalists. Spin Hall effect (SHE) [3, 4] and spin orbit 
coupling(SOC) are the most important candidates for many theoretical understanding in the realm of semiconductor 
spintronics. The spin Hall effect is the spin analogue of charge Hall effect with some differences as well. In spin 
Hall effect, external magnetic held is redundant for separating spin up and down spin electrons. The candidate 
responsible for the separation is SOC, which effectively generates a magnetic held in the rest frame of the electron 
and as a consequence we have spin current in this system. SOC, which is the relativistic coupling between orbital and 
spin degrees of freedom of electron can be obtained through the Foldy-Wouthuysen (FW) transformation [5] of the 
Dirac equation in presence of the external electric held. Alternatively, synthetic SOC can be generated via the strain 
parameter or with mechanical parameters like acceleration and rotation [6-8]. 

Besides, in semiconductor, the spin dynamics is inhuenced by the k.p perturbation theory [9]. The semiconductor 
band structure, close to band edges can be very well described by the k.p method. It is possible to explain the 
spin dynamics of semiconductor by taking into account the interband mixing via k.p perturbation theory. In [6], we 
have demonstrated that when band structure of semiconductor is considered, the free electron SOC parameter gets 
modihed by the Kane model parameters. The inclusion of this renormalized SOC parameter makes the theory of 
electron in semiconductor more accurate. 

Finite temperature effects are very important issues in different aspects of spin physics [10] . This gives the birth 
of spincaloritronics [II]. Very recently, the spin Hall conductivity is demonstrated in room temperature [12], which 
motivates us to study the thermally activated spin related issues. We are considering here the Fokker-Planck equation 
for analyzing the semi-classical motion of charge carriers. We have incorporated the additional constraints like damping 
force and also have included stochastic force arising due to the coupling of the system with a stochastic source of 
heat bath. In our formalism, the temperature correction is arising through the damping force. Besides, in presence of 
temperature the Kane model parameters are affected as well. This consequently affects the spin orbit coupling and 
electron ”g” factor. SOC is an important ingredient to have control over different physical parameters like spin current, 
conductivity. Berry curvature, spin relaxation time [13] e.t.c. In this paper, we theoretically have investigated the 
thermally driven spin current in a semiconductor on the basis of k.p perturbation theory from a generalized spin orbit 
Hamiltonian, which includes the stochastic force and arbitrary damping force. Here the Fokker-Planck or Kramer’s 
equation is employed to calculate the spin current. Furthermore, the effect of the crystal symmetry is also taken 
care of. At first, we have considered the temperature correction due to scattering mechanism through the damping 
constant 7. Secondly, the temperature dependence of the Kane model parameters is appraised. Our goal is to examine 
the expression of the thermally assisted spin current and spin conductivity in semi-conducting system. 

The organization of the paper is as follows: in section H we build our model Hamiltonian considering the k.p 
coupling between the Fg conduction band and Fg and Fy valance bands. In section HI, the semi-classical equation 
of motion is calculated applying the Kramer’s equation. We incorporate the effect of temperature through damping 
constant in this section and have computed the spin current. In IV, the renormalization of the Kane model parameters 
through temperature is taken care of, which modifies the SOC parameters as well. This renormalization of the SOC 
parameter alters the spin current in a different manner than that of the previous case. The conclusion is presented in 
section V. 
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II. THE MODEL HAMILTONIAN 


The Pauli-Schrodinger Hamiltonian with the effect of spin orbit coupling due to an external electric field can be 
written as [6-8] 


H = 


2m 


qU(r) + qX(7.{k X E) + 


( 1 ) 


where the first and second terms are the kinetic term with m as the free electron mass and the potential of the external 
electric field respectively. The third term is the spin orbit coupling term and forth term denotes the Zeeman term 
appearing as a consequence of external magnetic field. The free electron Hamiltonian in (1) modifies significantly 
when we consider the whole picture within a semiconductor, where one should incorporate the 8x8 Kane model [14] 
to include the effect of energy bands. The Hamiltonian for the 8x8 Kane model can be written as [6, 9] 
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and I 2 ,Ia are unit matrices of size 2 and 4 respectively. U = Ve(f) + Vc{r), is the total potential of the system which 
contains potential due to the external electric field Ve{f) and crystal potential Vc{r). Ec and Ey denote the energies 
at the conduction and valence band edges respectively. Ag is the spin orbit gap, P is the Kane momentum matrix 
element which couples s like conduction bands with p like valence bands. This Kane Momentum matrix element 
remains almost constant for group HI - V semiconductors, whereas Ag and Eq = E^ — Ey varies with materials. Here 
Eg denotes the energy gap between the conduction and valance band. The parameters P, Ag and Eg are known as 
the Kane model parameters. The Hamiltonian (3) can now be reduced to an effective Hamiltonian of the conduction 
band electron states [6, 9] as 
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(5) 

Now to find out the total Hamiltonian, we must add up the Hamiltonian (1) with (5). The total Hamiltonian for the 
electron in the conduction band edges can be written as [9] 


h2p ^ ^ So ^ 

Htot = 2 ^* 4" 4" 4- S\)a.{k x K) + (1 + —)/rtT .B, 


( 6 ) 


where 7* “ m 4- effective mass and E = —Vt4(r) is the effective total electric field and 

X = is the spin orbit coupling strength as considered in vacuum. Furthermore, the perturbation parameters JA 

and Sg are given by [9] 
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The 5X parameter is responsible for the renormalization of spin orbit coupling and the 5g term modifies the electron 
g factor considerably. It is possible to show that this extra term in the electron g factor can produce a shift in the 
ESR frequency. The Hamiltonian (6) can be rewritten neglecting the effect of Zeeman term as 
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where Ae// = A + <5A is the effective SOC term. The Hamiltonian in ( 8 ) is our system Hamiltonian, where the first 
term is the kinetic term, second term is the potential energy term and the third term denotes the SOC term. The 
renormalization of the mass and the SOC indicates that when we consider the electron within a semiconductor, we 
must take care of these Kane model parameters as well. 

The renormalized SOC parameter Ag// must influence the spin dynamics in of electron[ 6 ]. Our job is to find the 
spin current from equation ( 8 ). One can note that SOC is very important term in explaining the spin Hall effect. 
Here due to the interband mixing, the SOC term is changed. As a consequence the spin Hall current should modify 
as well. But how this SOC parameter is related to thermal corrections, is an important observation and we proceed 
to find this in section IV. But before that in section HI we want to find out the spin current without incorporating 
the thermal corrections of SOC. 


III. FOKKER-PLANCK EQUATION AND SPIN CURRENT 


Considering Hamiltonian ( 8 ), it is possible to calculate the semi-classical equations of motion by evaluating r and p 
via Heisenberg algebra. Before doing that, let us include the stochastic forces ({r, t), which appears as a consequences 
of other degrees of freedom as imperfection. One can also incorporate an arbitrary damping force K{r,p). Including 
all these forces, we can write the semi-classical equation of motion as 


• \ v 

r = TV [r, Htot] = ^ + eAe//(CT x VC/) 
in m* 


P = -^[p,Htot] =- eXeffV p.{axVU) - K{r,p) + C,{t). 


(9) 


As we are focusing in the linear response regime, it is suitable to consider C as distributed by Gaussian white noise with 
an arbitrary noise strength D. Thus the Fokker-Planck equation for the probability density P(r,p,t) can be obtained 
as [15] 
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The equilibrium solution can thus be written as 
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From definition, we can write k = 1 (p + m*Xeff{a x Vt/)^ , where 7 = 77 ^^- We have to find out the 

steady state solutions, K is the Boltzmann constant. Here we can write eqn. (10) as follows 
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where due to the presence of SOC, the momentum p is modified as 
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and A = — VC/— eAe//V p.{a xVU) — K{r,p). 

The time derivative of a as 

Including all these quantities the force equation has the following form 

(F) = = (A) + eXeffin.Via x VC/)). 

For a steady current, the average velocity of the carriers is fixed. The current in this case can be written as 
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FIG. 1: (Color online) Variation of the ratio of spin to charge conductivity with temperature for different semiconductors. 
Here a = . The unit in the x axis is Kelvin and in the y axis is A~^Kg~^ 


with the steady state constrain as 




( 17 ) 


Here n is the density of charge carriers and average is taken with the steady state solution of P. The steady state 
solution basically depends on the electric field and also on the coupling strength(Ae//). The expression for current 
for an isotropic crystal can be written as 


3 = ~, 

TO 7 

where the parameter /i is given as [16] 


E--{VVc} 
e 


eAe//ne^/i(l — a) 


TO* 7 ^ 


(ct X VU), 


d^Vc 

dr^dr^ 


= 




and a is a dimensionless parameter given as 


(VK,)=ae.E+i^(l-a)(JxE) 

Incorporating this in the expression (18) we have 

ne^(l — cx) 2Ae//ne^/x(l — cx) 

- - — 


J = 


m^"f 


TO*7^ 


(<7 X E) 


(18) 


(19) 


( 20 ) 


( 21 ) 


The first term corresponds to the charge current and the second term is the spin current. Thus the thermal effect is 
incorporated in the current via the 7 parameter. The expression of spin current can thus be written as 






( 22 ) 


where in the last step we have incorporated the value of 7 = This ensures that the spin current varies with 

temperature as T^. Similarly the charge current can be obtained as 
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Eqn. (22) and (23) suggests that the spin and charge conductivities is given by 
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This means that the spin to charge conductivity ratio is not constant rather depends on the temperature and 
increases with temperature. Also this ratio depends on the material chosen via the effective SOC parameter. Lastly 
one should also note that the above ratio also depends on the crystal symmetry parameter fj,. For anisotropic crystals 
this ratio should also change via the parameters associated with the crystal anisotropy [17]. FIG 1 shows the variation 
of the ratio of spin to charge conductivity with temperature for three different semiconductors. For the plot we have 
used the following table as 


EaieV) 

Ao(eK) 

P{eV°A) 

m*{me) 

GaAs = 1.519 

0.341 

10.493 

0.0665 

AlAs = 3.13 

0.300 

8.97 

0.150 

InP = 1.42 

0.110 

8.850 

0.0803 


As it is evident from the Fig. 1, that the dependence of the ratio on temperature is different for the three different 
semiconductors. This is one of the main results of this paper. As far as our knowledge goes, there is no experimental 
works which incorporate such stochastic forces to the system and investigate the effect of this on spin conductivity. 
Our work thus can initiate a possibility to experimentally verify the results stated here. 


IV. THERMALLY DRIVEN EFFECTS 

In the previous section, we have considered the effect of temperature on spin current, without incorporating the 
corrections due to thermal effects of SOC. In this section, our goal is to address the scenario of incorporating thermal 
consequences due to both stochastic force and modified Kane model parameters [14]. The temperature dependence 
of the band gap energy Eq can be given as [18] 


Eg = EoiO) - (26) 

where a,b are the fitting parameters or the varshni parameters [18], which are different for different semiconductors. 
The form of the energy band gap indicated in eqn. (26) is different from that of the expression in ref. [19], but they 
will lead to the same results, which is also mentioned in [19]. In addition to the temperature dependence of the gap 
parameters, the temperature dependence of the momentum matrix element P should also be considered. The Kane 
momentum matrix element P, varies with the lattice constant a as P « a(^' Here, the effect of phonon induced 
fluctuation of the interatomic spacing [19] is not included to avoid complexities in the calculations. The temperature 
dependence of the lattice constant can be written by the following relation [20, 21] 

a(r) = ai[l + at^(T- 300)], (27) 

where ath is the linear thermal expansion coefficient and the its value corresponds to the associated semiconductor. 
The values of Varshni’s parameters as well as of ath for two direct gap semiconductors can be given by [20, 22] 
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FIG. 2: (Color online) Variation of the ratio of spin to charge conductivity with temperature for GaP semiconductor, where 
a = The unit in the x axis is Kelvin and in the y axis is A~^Kg~^. 


is the effective mass. Although the spin orbit gap Ag may also be effected by the temperature, we are considering it 
as constant for simplicity. Thus incorporating eqn. (26) and (27) in eqn. (29), we have 
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Thus following the same approach as in the previous section, we have the expression of spin current as 
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. The spin current obviously depends on the temperature as well as the SOC strength A 
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where Ci = ^. The spin Hall conductivity varies with the temperature in a different manner. In the Previous case, 
it was a linear function of temperature. Figure 2 shows the variation of the dependence of the ratio of spin to charge 
conductivity with temperature. This is one of our main results in this paper. In order to show the variation of this 
spin Hall conductivity with temperature, we can write the above equation in the following form 

^ ^ (A + g(T)) Egg E^(0) 1 . W aT Ega^T ^ 

^ 1 'll I AbK^ IQbTK^ T^^SK^’ Ab'^K'^ 2bK 8K%‘^ J ’ ^ ’ 

writing the above equation only for the first order of Eq, a and b, we have 
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(A + C(r)) 
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rs 



Equ aT A 
AbK‘^ ■ 


(37) 


It is quite obvious from equation (37) that the spin conductivity depends on the temperature as well as the Varshni 
parameters a,b. One should note that, the temperature dependence is straightforwardly arising as a result of the 
thermal dependence of the Kane model parameters. We are not incorporating any temperature dependence forcefully. 
Interestingly, we are not considering any temperature gradient to show the thermal dependence of spin current. This 
is the beauty of our approach. 

We are now in a position to analyze the results of our paper. The expression of spin conductivity in eqn. (34) can 
be represented as 
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One can incorporate the second exponential to the constant term C and can rewrite the above expression as 
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where C = Cexpj^^y^jr^. At high temperature the exponential term dominates over the prefactor function of 
temperature. So we will achieve a exponential variation of the ratio. At low temperature the term C — — f, Ts is 

more dominant term. Also the result is strongly dependent of the chosen material parameters i.e on the Kane model 
parameters of the system. In these materials as the spin orbit gap parameter is small enough the key role is played by 
the energy parameter Eg and Kane momentum matrix parameter P, which have their variations with temperature. 


V. CONCLUSION 

In this paper, we have discussed the thermally driven effects on spin transport. Here, the k.p perturbation theory 
is employed to demonstrate the spin transport issues of an electron within a semiconductor. The charge carriers are 
influenced by the external electric field and a Gaussian white noise. The semi-classical force equation is calculated 
incorporating external forces due to damping and stochastic forces. The Kramer’s equation is adopted for evaluating 
the spin current. The temperature dependence appears due to the scattering mechanism induced in the damping 
factor 7 . 

It is well known that when we consider the electron transport within a semiconductor, we must take care of 
the band theory of the semiconductors. The inclusion of the band theory, enables us to include the Kane model 
parameters in the theory. In presence of temperature, the Kane model parameters are perturbed. This consequently 
affects the SOC parameters. The temperature dependence of the SOC parameter gives a thermally dependent spin 
current. The thermal dependence of spin Hall conductivity is also discussed in this paper. This shows a different 
approach of encountering the thermal dependence scenario in the spin transport regime. 
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